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ABSTRACT 

Be stars are a class of rapidly rotating B stars with circumstellar disks that cause Balmer and 
other line emission. There are three possible reasons for the rapid rotation of Be stars: they may 
have been born as rapid rotators, spun up by binary mass transfer, or spun up during the main- 
sequence (MS) evolution of B stars. To test the various formation scenarios, we have conducted 
a photometric survey of 55 open clusters in the southern sky. Of these, five clusters are probably 
not physically associated groups and our results for two other clusters are not reliable, but we 
identify 52 definite Be stars and an additional 129 Be candidates in the remaining clusters. We use 
our results to examine the age and evolutionary dependence of the Be phenomenon. We find an 
overall increase in the fraction of Be stars with age until 100 Myr, and Be stars are most common 
among the brightest, most massive B-type stars above the zero-age MS (ZAMS). We show that 
a spin-up phase at the terminal-age MS (TAMS) cannot produce the observed distribution of Be 
stars, but up to 73% of the Be stars detected may have been spun-up by binary mass transfer. 
Most of the remaining Be stars were likely rapid rotators at birth. 

Previous studies have suggested that low metallicity and high cluster density may also favor Be 
star formation. Our results indicate a possible increase in the fraction of Be stars with increasing 
cluster distance from the Galactic center (in environments of decreasing metallicity). However, 
the trend is not significant and could be ruled out due to the intrinsic scatter in our data. We 
also find no relationship between the fraction of Be stars and cluster density. 

Subject headings: stars: emission-line, Be — open clusters and associations: individual (Basel 1, Bochum 
13, Collinder 272, Haffner 16, Hogg 16, Hogg 22, IC 2395, IC 2581, IC 2944, NGC 2343, NGC 2362, 
NGC 2367, NGC 2383, NGC 2384, NGC 2414, NGC 2421, NGC 2439, NGC 2483, NGC 2489, NGC 
2571, NGC 2659, NGC 3293, NGC 3766, NGC 4103, NGC 4755, NGC 5281, NGC 5593, NGC 6178, 
NGC 6193, NGC 6200, NGC 6204, NGC 6231, NGC 6249, NGC 6250, NGC 6268, NGC 6322, NGC 
6425, NGC 6530, NGC 6531, NGC 6604, NGC 6613, NGC 6664, Ruprecht 79, Ruprecht 119, Ruprecht 
127, Ruprecht 140, Stock 13, Stock 14, Trumpler 7, Trumpler 18, Trumpler 20, Trumpler 27, Trumpler 
28, Trumpler 34, vdB-Hagen 217) 

1. Introduction 

1 Visiting Astronomer, Cerro Tololo Inter-American Ob- 
servatory. CTIO is operated by AURA, Inc. under contract Be stars are well known to be a class of rapidly 
to the National Science Foundation. rotating stars (Slettebak 1949, 1966; Slettebak et 

2 NSF Astronomy and Astrophysics Postdoctoral Fellow gj 1992). Slettebak (1966) showed that the ob- 
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served equatorial rotational velocities, v rot , of Be 
stars are 70%— 80% of the critical breakup veloc- 
ity, Vcrit, although Townsend, Owocki, & Howarth 
(2004) and Fremat et al. (2005) showed that ob- 
served Be star rotational velocities may be system- 
atically underestimated due to the effects of gravi- 
tational darkening on spectral lines. Therefore Be 
stars may rotate much closer to their critical veloc- 
ities than the available observations suggest. This 
rapid rotation may be combined with weaker pro- 
cesses, such as non-radial pulsations or magnetic 
fields (Porter & Rivinius 2003), to move material 
from the stellar surface into a disk (Townsend et 
al. 2004). 

There are three possible reasons for the rapid 
rotation of Be stars: they may have been born as 
rapid rotators, spun up by binary mass transfer, or 
spun up during the MS evolution of B stars. Many 
authors (Mermilliod 1982; Grebcl 1997; Zorec & 
Briot 1997; Fabregat & Torrejon 2000) have found 
the highest fractions of classical Be stars among 
spectral types B0— B2, appearing most often be- 
tween the ages 13 — 25 Myr (Fabregat & Torrejon 
2000). Fabregat & Torrejon use these results to 
argue that the Be phenomenon is due to an evo- 
lutionary spin-up among B-type stars towards the 
end of the MS lifetime (Meynet & Maeder 2000; 
Keller et al. 2001). On the other hand, Abt (1979) 
found no dependence on the Be frequency with 
age, and Mermilliod (1982) and Slettebak (1985) 
noted that Be stars are found throughout the en- 
tire MS band, from the ZAMS through the TAMS, 
indicating that Be stars can occur at any age. 
Therefore it is also possible that the evolution of 
close binary systems may lead to mass transfer 
onto B stars, thereby increasing their angular mo- 
mentum and inducing disks (Pols et al. 1991). Fi- 
nally, Zorec & Briot (1997) find constant fractions 
of Be stars with luminosity class, and they argue 
that the Be phenomenon is not associated with a 
particular evolutionary phase. Rather, it is more 
likely due to rapid rotation in B stars at the time 
of their formation. Because these various stud- 
ies disagree about the relationship between the Be 
frequency and evolution, further investigations are 
necessary to resolve the issue. 

Therefore we present here an extensive pho- 
tometric survey of open clusters to study their 
Be star populations. Our detection technique is 
demonstrated for the cluster NGC 3766 in Mc- 



Swain & Gies (2005), hereafter referred to as Pa- 
per 1. In this paper we present photometric data 
for 55 open clusters, and we identify a total of 52 
definite Be stars and 129 Be star candidates in 
the 48 reliable open clusters. We use our sample 
of Be stars to conduct a statistical analysis of the 
Be phenomenon and explore the reasons for their 
rapid rotation. 

2. Observations and Data Analysis 

We selected the target open clusters for this 
study using the WEBDA database 3 . We chose 
clusters with ages between 1-32 Myr, V — My < 
13, and angular diameters between 5-20". How- 
ever, we also included some older clusters to im- 
prove our statistics, and many of the ages and 
distance moduli from WEBDA were later revised 
using other references. We restricted our targets 
to clusters in the southern sky; all have declina- 
tion < 0°. To eliminate clusters with obvious as- 
sociated nebulosity, we used the Space Telescope 
Science Institute (STScI) Digitized Sky Survey 4 
to view a red image of each cluster. While the 
requirement of little to no nebulosity reduced the 
difficulty of the photometry, it did eliminate many 
very young clusters from our sample. 

We made photometric observations of the clus- 
ters over 14 nights between 2002 March 27 - April 
2 and 2002 June 19 - 25 with the CTIO 0.9m 
telescope and SITe 2048 CCD. The images were 
binned using a CCD summing factor of 2 pixels 
x 2 pixels due to the slow readout time of the 
chip. Without binning, the chip has a plate scale 
of 0'.'401 pixel -1 , but with binning the plate scale 
increased by a factor of 2. 

Most clusters were observed at least twice in 
each of the Stromgren b and y filters as well as 
a narrow band Ha filter, and the exposure times 
usually varied from 5-120 s duration depending 
on the brightness of the cluster. During the 2002 
March run, we also observed five standard stars 
from the list of Cousins (1987): HD 79039, HD 
80484, HD 104664, HD 105498, and HD 128726. 
We observed the standard stars HD 105498, HD 

3 The WEBDA database is maintained by J.- 
C. Mermilliod and is available online at ob- 
swww.unige.ch/webda/navigation.html. 

4 The STScI Digitized Sky Survey is available online at 
http://stdatu.stsci.edu/cgi-bin/dss_form. 
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128726, HD 156623, HD 157795, HD 167321, and 
HD 216743 from the lists of Cousins (1987) and 
Clausen et al. (1997) during the 2002 June nights. 
Unless the weather prohibited observations, each 
standard star was observed in each band at a min- 
imum of three different airmasscs. All of the im- 
ages were processed in IRAF 5 using nightly bias 
and dome fiat frames. 

The photometry of each standard star and clus- 
ter was measured and calibrated according to the 
the methods discussed in Paper 1. We measured 
the errors in the magnitudes by combining the es- 
timated instrumental errors, the standard devia- 
tion of the aperture correction, and the errors in 
the transformation coefficients, and these quanti- 
ties vary for each cluster due to nightly variations 
in atmospheric conditions and the exposure time 
used. From our images, we also performed astrom- 
etry of the clusters according to the technique de- 
scribed in Paper 1, and the errors are typically less 
than O'.'IO in both acos<5 and 5. However, several 
clusters did not fit into a single field of view, and 
they were observed in three or four fields. These 
clusters are Collinder 272, Hogg 16, IC 2395, NGC 
2343, NGC 6425, NGC 6531, Stock 13, and Trum- 
pler 18. The astrometry errors of these clusters are 
larger, but usually less than 0750 in both a cos 5 
and 6. 

Our photometry and astrometry of each star 
in the 55 clusters are listed in Table 1, available 
online (see Note to Table 1). Although we often 
exclude stars near the edge of the CCD or near 
bad columns, the photometry of each open cluster 
is usually complete enough to include all of its 
(unsaturated) OB stars. It is complete to y = 13 
for the brightest clusters, but we obtained deeper 
images and measured magnitudes up to y = 17 in 
some cases. 

From the calibrated photometry, we created 
(b — y, y) color-magnitude and (6 — y,y — Ha) color- 
color diagrams like those described in Paper 1 for 
NGC 3766. Because the MS of the (b - y, y) dia- 
gram is nearly vertical in the region of massive OB 
stars, it was very difficult to constrain the distance 
modulus and the age of each cluster. Therefore, 

5 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



to determine the cluster parameters (E(B — V), 
Vq — My, and age), we searched the available lit- 
erature for the various values for each cluster. We 
used these values to plot the theoretical isochrones 
from Lcjcunc & Schaerer (2001) and theoretical 
color-color curve using colors derived from Kurucz 
(1979) model spectra. This technique is described 
for NGC 3766 in Paper 1. The Lejeune & Schaerer 
(2001) isochrone models do not provide Stromgren 
magnitudes, so we transformed the Johnson V and 
B magnitudes to the Stromgren system using the 
transformations 

y = V - 0.038(5 - V) (1) 

(Cousins & Caldwell 1985) and 

B- V = 1. 584(6 -j/) + 0.681toi -0.116 (2) 

(Turner 1990), as described in Paper 1. We also 
transformed E(B — V) to E(b — y) using the rela- 
tionship 

E(b - y) = 0.745 E(B - V) (3) 

from Fitzpatrick (1999), assuming R = A(V) /E(B— 
V) = 3.1. In very few CRS6S, cL cluster did not have 
published parameters available, or the published 
values did not agree well with our observations, so 
we fit our data to the theoretical curves using a 
grid of values for E(b — y), Vo — My, and log age. 
For all clusters, we retained the parameters that 
resulted in the best fit between our data and the 
theoretical curves, and the resulting parameters 
are listed in Table 2. Our technique to select defi- 
nite and possible Be stars from the normal B stars 
in our images according to their excess brightness 
in the narrow band Ha filter is described in Paper 
1, and Table 2 also includes the numbers of Be 
and B stars identified in each cluster. 

3. Results from Photometric Study 

Before we consider the results from this photo- 
metric survey, it is important to review the errors 
inherent in this study. Saturated stars affect the 
number of very bright B and Be stars that we de- 
tect, but this problem is ameliorated by the avail- 
ability of data on most bright cluster members in 
the literature. However, we do not include sat- 
urated stars in the number of Be and B stars in 
Table 2. A second problem is that our photometry 
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technique only detects those Be stars that had ac- 
tive Ha emission at the time of our observations, 
so we detect only the lower limit of the true Be 
population. We showed in Paper 1 that the true 
fraction of Be stars is probably between 1.2 — 2.3 
times the number detected. Assuming that the ra- 
tio of active to inactive Be stars is similar among 
different clusters, the neglect of inactive Be stars 
will have no impact on our statistical comparison 
of the Be star frequency between clusters. The 
largest source of errors in our study is due to the 
difficulty detecting weakly emitting Be stars since 
the errors in their photometry are comparable to 
the strength of their emission. There is also a 
possibility that unreddened foreground stars or, 
less often, background supcrgiants, contaminate 
the MS and appear as weak emitters in the color- 
color diagram. Field stars cause special problems 
for the five targets that arc not trucly bound clus- 
ters, and the ages, reddening, and distance moduli 
for these groups are not representative of the en- 
tire stellar population. Finally, one of our require- 
ments for identifying B-type stars (see Paper 1) is 
that b — y < E(b — y), which neglects the inherent 
width of the MS, and we may omit some B and 
Be stars redward of this cutoff. This is especially 
a problem at the faint end of the B star MS. 

Our photometry of six clusters, excluding the 
less reliable clusters Bochum 13 and Trumpler 27, 
did not cover the entire range of B star magnitudes 
so the total number of B stars in these clusters is 
underestimated in Table 2. We show below that 
the number of Be stars omitted is negligible, in 
agreement with the results of Keller et al. (2001). 
However, we account for the missing late B-type 
stars by computing the mass function, £(logM). 
The Lejeune & Schaerer (2001) isochrones assign 
a mass to each point along the curve, and we in- 
terpolated to convert the observed y magnitudes 
of the B stars to M. For the allowed mass range 
of B-type stars, we binned the stars according to 
logM and counted the unnormalized number, N, 
in each bin. The slope of the linear fit between 
log TV and logM provides F, where 

r = dlog<e(logM)/dlogM. (4) 

In the fit, each value of N was weighted according 
to TV 1 / 2 , which decreases the dependence on the 
less common high mass stars, and the incomplete 
lower mass bins were excluded entirely. The aver- 



age log M in each bin was used to compute the to- 
tal number of B-type stars expected, N tot . We cal- 
culated T and N to t using 6—15 log M bins, equally 
spaced along the MS curve in the color-magnitude 
diagram. Table 3 provides T, N toU N missing , and 
the observed mass range of B-type stars, AM, for 
the six clusters analyzed. This mass range varies 
depending on the cluster age, and the upper mass 
limit is comparable to the observed maximum B 
star mass in each cluster. The minimum mass of 
B stars is 2.3 M Q from the Lejeune & Schaerer 
(2001) models. We found a slight dependence be- 
tween r and the number of bins used, and the 
errors for T and N to t reflect the la variance. Our 
values for T are consistent with the values for mas- 
sive stars in open clusters found by Massey et al. 
(1995). 

In this study, we find a total of 52 definite Be 
stars and an additional 129 possible Be stars in 48 
clusters. (Hereafter we ignore the regions that are 
not true open clusters as well as Bochum 13 and 
Trumpler 27, discussed in Appendix A.) We also 
detect 2210 B-type stars, including the Be stars, 
and we estimate from the clusters' mass distribu- 
tions that 328 B-type stars are missing from our 
sample. Considering only the definite Be stars in 
these clusters, the mean percentage of Be stars is 
only 2.0%. Including the possible Be stars raises 
the mean percentage to 7.1%. This is much lower 
than expected; Abt (1987) found that 18% of field 
B0 to B7 III-V stars in the Bright Star Catalogue 
(BSC; Hofflcit & Warren 1991) were classified as 
Be stars in the literature. However, our study ac- 
curately accounts for the total number of B stars 
within the cluster volume, while Abt's magnitude- 
limited survey is biased towards early B stars that 
are probably more rich in Be stars (e.g. Mermil- 
liod 1982). Furthermore, many of the B stars in 
the BSC are members of the Gould Belt, which 
has 7.4 < log age < 7.8 (Perrot & Grenier 2003) 
and coincides with a maximum in the Be star frac- 
tion (sec Figure 5 below). This may explain the 
unusual richness of Be stars in the BSC. 

Using our extensive sample of Be stars, we in- 
vestigate several factors that may influence the 
population of Be stars in open clusters, begin- 
ning with metallicity. Using clusters within the 
same age interval in the interior and exterior of 
the Galaxy, the Large Magellanic Cloud (LMC), 
and the Small Magellanic Cloud (SMC), Maeder, 
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Grebel, & Mermilliod (1999) found that the ra- 
tio of Be stars decreases sharply with increasing 
metallicity. However, Keller, Wood, & Bessell 
(1999) found no difference in the fraction of Be 
stars in the LMC and SMC, although there is 
a factor of two difference in their metallicities. 
To look for this influence in our data, we inves- 
tigated the distribution of the target clusters in 
the Galactic plane. We calculated the distance of 
each cluster from the Galactic center, d cen t, us- 
ing R — 8000 pc as the Sun's distance from the 
Galactic center and assuming that each cluster lies 
within the plane of the Galaxy (valid since all have 
| b |< 6°). The Galactic longitude, £, latitude, b, 
distance from the Sun, d Q , and d cen t for each clus- 
ter are listed in Table 4. The spatial distribution 
of clusters with respect to the Sun and the Galac- 
tic center is illustrated in Figure 1. 

The numbers of definite and possible Be star 
detections from our photometry arc included in 
Table 2, and we plot the percentage of Be stars 
as a function of distance from the Galactic cen- 
ter in Figure 2. The diamonds are the lower lim- 
its of the Be star fraction, determined from the 
number of definite detections, and the upper limit 
assumes that all of the possible candidates are in 
fact Be stars. Note that this upper limit does not 
represent the formal error in our measurement, 
which can be significant in clusters with few B- 
type stars (especially NGC 6268, in which we find 
that one of its two B stars is a Be star). For the 
35 clusters with d cen t < Ro, the mean Be star 
population is 1.8%— 6.5%. The outer 13 clusters 
have a mean population of 3.3%-10.5%. Using the 
Galactic metallicity gradient AZ/kpc = —0.0019 
to -0.0029 (Maeder et al. 1999), over the 4.3 
kpc range in distance from the Galactic center, 
AZ = .0082 - .0125. Interpolating from Maeder 
ct al., the fraction of Be stars in our clusters should 
increase by a factor of « 2 over this range in metal- 
licity. The mean Be star populations do increase 
slightly over this distance, but Figure 2 does not 
provide compelling evidence supporting a depen- 
dence on the Galactic Be star fraction with metal- 
licity. 

The number of Be stars may also be related to 
the stellar density in the cluster. Strom, Wolff, 
& Dror (2005) suggest that the average rotational 
velocity is greater for higher density clusters, pre- 
sumably allowing more Be stars to form. They 



show that the average rotational velocity of un- 
evolved B stars in the dense clusters h and \ Persei 
is more than 2 x that of field stars of comparable 
ages. To investigate the role of the stellar den- 
sity in each cluster, we calculated the mean B star 
density of the clusters. We determined the clus- 
ter center using the mean positions of the B-type 
stars, and the standard deviation of the B star 
positions relative to the cluster center defines the 
radius of the assumed spherical volume. However, 
the cluster NGC 6268 has only two B stars, so we 
used the mean distance from the cluster center to 
define the radius in this case. We plot the per- 
centage of Be stars against the number of B stars 
per cubic parsec in Figure 3. Our data show no 
correlation between the two. 

We also consider the relationship between the 
number of Be stars in open clusters and the clus- 
ter age, a possible signature of evolutionary ef- 
fects involved in the origin of Be star disks. Past 
studies have disagreed on the relationship between 
the ages and numbers of Be stars. Fabregat & 
Torrejon (2000) found that classical Be stars are 
most abundant between 7.1 < log age < 7.4. On 
the other hand, Mermilliod (1982) and Slettebak 
(1985) found a weaker dependence on age. In- 
stead, they noted that Be stars occupy the entire 
MS band, from the ZAMS through the TAMS, 
indicating that Be stars occur at any evolution- 
ary stage. Mermilliod (1982) found that the fre- 
quency of Be stars depends instead on spectral 
type; the abundance peaks for types B1-B2 and 
B7-B8. However, Zorec & Briot (1997) also found 
no distinct peak the Be star frequency with lu- 
minosity class, and they argue that the Be phe- 
nomenon is more likely due to the characteristics 
of B stars at the time of their formation. This 
would imply a higher fraction of Be stars among 
young clusters. 

We present a plot of the Be star percentage as 
a function of cluster age in Figure 4. While we 
did not compute the formal error in the ages, the 
agreement among different ages provided in the 
literature and the good match between the the- 
oretical isochrones and our color-magnitude dia- 
grams suggest that the cluster ages are accurate 
to ±0.2 dex in most cases. Because no clear re- 
lationship between the percentage of Be stars and 
the cluster age emerges from Figure 4, we plot the 
average Be frequency in four age bins in Figure 
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5. The bins correspond to the ages log age < 7.0, 
7.0 < log age < 7.4, 7.4 < log age < 8.0, and 
log age > 8.0. Our definite Be star frequency and 
the Poisson errors are plotted using solid lines, 
and the total of definite and possible Be stars, 
also with the Poisson errors, are represented with 
dashed lines. The plotted ages represent the mean 
age observed within each bin, and the age er- 
ror bars represent the maximum and minimum 
ages from our sample. We are surprised to see 
that the fraction of Be stars is generally highest 
among clusters with 7.4 < log age < 8.0, unlike 
previous studies that found a peak among clus- 
ters with 7.0 < log age < 7.4 (Fabregat & Tor- 
rejon 2000). Figure 5 reveals that the popula- 
tion of Be stars generally increases with time until 
7.4 < log age < 8.0, then declines among older 
clusters. 

Finally, we plot the 52 definite Be stars in 
Figure 6 using their intrinsic colors and absolute 
magnitudes. For comparison, the normal B-type 
stars are plotted as contours overlying the absolute 
color-magnitude diagram. To derive the absolute 
magnitudes, we used the relation 

A(y) = 3.034 E(B - V) (5) 

interpolated from Fitzpatrick (1999). The hori- 
zontal dashed lines in Figure 6 divide the Be pop- 
ulation into four bins: evolved (M y < yBmin at 
the ZAMS), early B spectral types, mid B stars, 
and late B stars. The isochrone for 10 5 yr repre- 
sents a very young age, close to the ZAMS, and 
the isochroncs for a distribution of stars with ages 
of 10 7 and 10 8 yr are also plotted. Although the 
Be stars are generally a redder population than 
the B-type stars (the two most significant outliers 
to the lower left of the ZAMS are No. 12 in NGC 
6531 and No. 355 in Trumpler 34; see discussion 
in Appendix A), it is unlikely that the Be stars 
in a given cluster are a more evolved population 
than their neighboring B stars, and we believe that 
other effects are largely responsible for their posi- 
tion relative to the ZAMS. The circumstellar disk 
can impart additional reddening, rotationally in- 
duced gravitational darkening may cause them to 
appear cooler than a non-rotating star of the same 
mass, or light may be preferentially scattered by 
the disks and cause the stars to appear brighter 
and more evolved (Slettebak 1985; Zorec & Briot 
1997). 



Although Figure 6 reveals that the definite Be 
stars are rather uniformly distributed in M y , the 
relative fraction of Be and B stars decreases with 
later spectral tyes. Be stars comprise a respec- 
tive 7.4%, 4.7%, 2.1%, and 0.9% of the evolved, 
early, mid, and late type B stars in our sample. 
A histogram plot illustrating the relative popula- 
tions for each bin is shown in Figure 7. These 
results indicate that the relative fraction of Be 
stars attains a maximum among the evolved stars 
that are approaching the TAMS. (Note that these 
evolved stars do not correspond to the oldest age 
bin defined above, where a decline in the Be pop- 
ulation is observed.) Fabregat & Torrejon (2000) 
also found that the Be phenomenon occurs dur- 
ing the second half of the MS, which they argue is 
the result of an evolutionary spin-up as the stars 
approach the TAMS. 

4. Conclusions 

Our results indicate that the Be phenomenon 
is not strongly dependent on mctallicity or clus- 
ter density, but it is clearly influenced by the evo- 
lutionary state of B stars. We observe a higher 
fraction of evolved Be stars, but if Be stars were 
born as rapid rotators as believed by Zorec & Briot 
(1997), the opposite trend would be observed; 
more Be stars would be observed in unevolved 
stars, in the youngest clusters, since the rotation 
rate decreases during the MS lifetime (Meynet & 
Maeder 2000). Also, Be stars born as rapid ro- 
tators would have a much different distribution 
along the MS. A rotationally distorted star with 
mass Mb at the critical breakup velocity will have 
an equatorial radius 

^ = 1 K - 914+ (^F) i?S ^ 

(Fremat et al. 2005), where Rb is the radius of a 
non-rotating B star. The resulting critical breakup 
velocity is 

vent = ^ (^^) ^ km s" 1 . (7) 

Table 5 compares this v crit with the mean pro- 
jected rotational velocity, (v sinz), of B stars from 
Abt, Levato, & Grosso (2002). Masses and radii 
are from Harmanec (1988). Later B-type stars 
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have a mean rotational velocity closer to their crit- 
ical value, so that if these velocities reflect their 
birth spins, then we should expect to see a higher 
fraction of Be stars among late B stars. We show 
in Figure 7 that such a distribution is not ob- 
served, hence we doubt that many Be stars are 
born as fast rotators. 

Meynet & Maeder (2000) predicted that stars 
with M > 12M experience a short spin-up phase 
at the end of their MS lifetimes, and this spin-up 
phase could be responsible for the high fraction 
of evolved Be stars observed in this study. We 
do not expect an age dependence with this scene- 
rio; as long as a cluster has TAMS B stars, Be 
stars would be expected. However, the duration 
of this spin-up phase is expected to be short, only 
about 1% of the MS lifetime, so we should observe 
a comparable Be star fraction if the TAMS spin- 
up phase is the origin of the rapid rotation. While 
lower mass B stars are not expected to experience 
such a spin-up, Meynet & Maeder (2000) predict 
that these stars may remain rapid rotators and in 
fact increase their fraction of the breakup veloc- 
ity as they evolve along the MS. Once again, this 
would support a larger fraction of late Be stars, 
which is not observed. 

We find it more likely that Be stars are spun- 
up by mass transfer in close binary systems. Al- 
though there is a perceived deficiency of short pe- 
riod Be star binaries (e.g. Abt & Levy 1978), re- 
cent work has led to the discovery of 55 Be star bi- 
naries with periods < 100 d (Gies 2000; Raguzova 
& Popov 2005), and Be stars have been observed in 
nearly all stages of the proposed binary evolution 
scenario. The initially more massive star expands 
to fill its critical Roche lobe during the H shell 
burning phase, triggering Case B mass transfer 
and causing the spin-up of the mass gainer (Gies 
2000). The massive Algol binary RY Per is an ex- 
ample of such a spin-up process in action; the F7: 
II-III donor star is engaged in active mass transfer 
via Roche lobe overflow onto the B4: V compan- 
ion, and the mass gainer has been spun up to more 
than seven times faster than the synchronous rate 
(Barai et al. 2004). Commonly, the mass ratio be- 
comes inverted during the mass transfer and the 
orbital separation a increases, leading to a dra- 
matic increase in the timescale for tidal synchro- 
nization (~ a 17 / 2 ; Hilditch 2001). This timescale 
becomes much longer than the H shell burning 



phase, so the mass gainer will retain its rapid ro- 
tation long after the mass transfer is complete. 
At the conclusion of the mass transfer phase, the 
donor star is left as a stripped down, helium star 
remnant - a small but hot subdwarf (Gies 2000). <fi 
Per is the only confirmed Be+sdO system (Thaller 
et al. 1995; Gies et al. 1998), but three additional 
candidates have been identified: 59 Cyg (Maintz et 
al. 2004), HR 2142 (Waters, Cote, & Pols 1991), 
and FY CMa (Rivinius et al. 2004). Depending 
on its mass, the sdO star will eventually become 
a white dwarf (WD) or explode as a supernova 
to form a Be + neutron star (NS) binary or an 
unbound system (Gies 2000). To date, no Be + 
WD systems are known, but 130 Bc/X-ray bina- 
ries (presumably with NS companions) have now 
been identified (Raguzova & Popov 2005). 

Pols et al. (1991) and Van Bever & Vanbeveren 
(1997) predicted that the fraction of mass gain- 
ers is highest among B0-B3 stars, decreasing with 
later spectral types, consistent with the distribu- 
tion of Be stars that we observe. If the mass trans- 
fer is conservative (as expected for systems with a 
high initial mass ratio, q > 0.6), Pols et al. (1991) 
showed that the mass gainers would be observed as 
rejuvenated, bright MS stars, including blue strag- 
glers. However, Be stars formed this way may not 
appear near the ZAMS due to disk reddening and 
gravitational darkening. On the other hand, Pols 
et al. showed that there is a minimum q (between 
0.3 — 0.5, depending on their model) for the pro- 
duction of Be stars from close binary evolution 
with non-conservative mass transfer. Therefore, 
Be stars formed by mass transfer would be found 
preferentially in the upper part of the MS. 

Previous studies have estimated the percent- 
age of the Be stars that are products of close bi- 
nary evolution. Pols et al. (1991) predicted that 
1.3%— 3.9% of B-typc stars should be observed 
as rapid rotators due to angular momentum re- 
ceived during case B mass transfer, and Van Bever 
& Vanbeveren (1997) predicted a similar number 
(0.9%-4.8%). While the predicted total number 
of mass gainers is comparable to the number of Be 
stars that we observe, not all of these mass gainers 
are expected to become Be stars. Pols et al. (1991) 
showed that the predicted number of spun-up B 
stars in the BSC is 40%— 60% of the observed num- 
ber of Be stars in that sample. However, because 
not all rapidly rotating B-type stars exhibit the Be 
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phenomenon, Van Bever & Vanbeveren (1997) ar- 
gued that Be stars represent about half of the total 
number of rapidly rotating B stars, so the fraction 
of Be stars in the BSC formed by mass transfer is 
closer to 20%. We show above that the fraction 
of Be stars in the BSC is not representative of the 
fraction in open clusters, so the estimate of Van 
Bever & Vanbeveren (1997) cannot be applied to 
our sample, and there is no spectroscopic survey 
in the literature that provides relative numbers of 
Be stars and normal, rapidly rotating B stars. 

In lieu of an accurate theoretical prediction, we 
can estimate the fraction of mass gainer Be stars 
from our results. We plot the observed M y of our 
definite Be stars versus their age in Figure 8. Us- 
ing the isochrone models of Lejeune & Schaerer 
(2001), we determined the mass of a primary at the 
TAMS and its M y (dashed curve) over a range of 
cluster ages. For a close binary system with q « 1 
and the primary at the TAMS, the mass transfer 
will be conservative (Pols et al. 1991) and the re- 
juvenated secondary will appear brighter than the 
TAMS My. Assuming a minimum q — 0.3 for the 
production of Be stars (Pols et al. 1991), we also 
identified the M y of the mass gainer, shown as a 
dotted curve. (Note that for log age > 7.5, the lim- 
iting My for the secondary is greater than ysmax-) 
Figure 8 reveals that only one of our definite Be 
stars is found at the TAMS and thus may have ac- 
quired its rapid rotation during a spin-up phase at 
this stage of its evolution. In contrast, we find that 
38 Be stars, or 73%, fall in the mass range associ- 
ated with binary mass ratios of 0.3 < q < 1 and 
are good candidates for spin-up by binary mass 
transfer. (Note that two Be stars in NGC 2439 
have nearly the same M y « —2.1, so these points 
overlap in Figure 8.) There are 14 Be stars below 
the binary mass ratio line of q < 0.3, and some of 
these were probaby born as rapid rotators. How- 
ever, four of these (members of the clusters NGC 
6231, NGC 6322, and Stock 13) are likely Her- 
big Be stars, and their disks may be due either 
to rapid rotation or the remnants of star forma- 
tion. We also caution that the M y relations in 
Figure 8 neglect gravitational darkening and disk 
reddening in Be stars, so our comparison of the 
predictions of the three Be formation scenarios is 
not exact. Nevertheless, our results clearly empha- 
size the relative dominance of the binary spin-up 
process, while neither a TAMS spin-up phase nor 



rapid rotation at the time of formation appear to 
play a significant role in the formation of Be stars. 
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A. Notes on Individual Clusters 

While the results from each cluster are summarized by Table 2, further comments regarding each cluster are 
presented below. The definite and possible Be stars in each cluster are identified using their assigned number 
from Table 1. (For cross referencing purposes, note that Table 1 lists other identifiers and WEBDA numbers 
for these stars whenever possible.) Particular sources of error, including saturated stars and contaminating 
foreground stars, are also discussed. Unless otherwise noted, spectral types and other information about 
individual stars are from the SIMBAD database. 

The (b — y, y) color-magnitude and (b — y, y — Ha) color-color diagrams for each cluster are presented in 
Figures 9 — 63. In each color-magnitude diagram, the closest matching theoretical isochrone, assuming solar 
metallicity, from Lejeune & Schaerer (2001) is plotted as a solid line. The limiting magnitudes of B-type 
stars, ysmin and yBmax, are indicated as dashed lines when they are within the plotting range. In each 
color-color diagram, the theoretical color-color curve, derived from Kurucz (1979) model spectra, is shown 
as a solid line. Be stars are identified by their excess Ha brightness that shifts their y — Ha color above 
this reference line (Paper 1). We also made a parabolic fit of the color indices derived from dereddencd 
spectra of 161 stars from the atlas of Jacoby, Hunter, & Christian (1984), shown as a dashed curve, and 
unreddened foreground stars will probably be found close to this curve in the observed color-color diagrams. 
In both diagrams, the massive O- and B-type stars all lie to the left of the vertical dashed line that represents 
E(b — y). The definite and possible Be stars are represented with large diamonds in each plot, while all other 
stars are represented with small diamonds. In a few cases, potential foreground stars or white dwarfs that 
lie to the left of the reddening line have been excluded from the Be sample, and these stars are marked with 
an X. 

• Basel 1: Only one star, GSC 05126 — 02385, is saturated in our images of Basel 1, and we detect no 
Be stars in this cluster. 

• Bochum 13: Our images of Bochum 13 contain two bright, saturated stars: HD 156134 and HD 
156154. The former is an 0+ star, and the latter is a B0 lab star. We identify Nos. 16 and 30 as 
definite Be stars, and No. 22 is a possible Be star candidate. 

Aside from the two bright stars that are saturated, Bochum 13 is a sparse cluster that contains mostly 
fainter members. The MS is not well defined, and we find it difficult to determine accurate estimates 
of E(b — y), Vo — My, and log age for this cluster. Therefore we consider our results less reliable, and 
we ignore this cluster in our statistics. 

• Collinder 272: There are no saturated stars in our images of the cluster Collinder 272. The three Be 
stars in this cluster are Nos. 508, 565, and 767. We also find five possible Be stars: Nos. 345, 513, 560, 
562, and 774. 

• Haffner 16: In our images of Haffner 16, the cluster appears to have a small diameter of approximately 
4', and many field stars are visible in the region. However, to create the color-magnitude and color-color 
diagrams of the cluster, we used the stars from the entire 13.'5 x 13f 5 field of view of the CCD. From 
our color-magnitude diagram shown in Figure 12, the region of sky appears to contain two distinct 
stellar populations. Table 4 and Figure 1 reveal that Haffner 16 is located within the Perseus arm of 
the Galaxy and is viewed through the local spur, causing significant contamination from foreground 
stars along the line of sight. 

Four stars in this cluster show Ha emission along with a color b — y < 0.123. One of the stars, No. 
756, is within the dense cluster region and we classify it as a possible Be star. The others have a y 
magnitude fainter than the B star limit for the cluster, so they are likely foreground stars. There are 
three stars in our images of Haffner 16 that are saturated, but none are identified in the SIMBAD 
database. 
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Hogg 16: Two stars are saturated in Hogg 16: HD 116887, a KO III star, and HD 116875, a B8 V 
star. We find one possible Be star in Hogg 16, No. 195. The SIMBAD database also identifies No. 199 
as a B2e star, and while we do observe Ha emission from this star, its b — y color is too red to classify 
it a B-type star in the cluster. 

The cluster Hogg 16 lies about 16' west of Collinder 272. We find slightly different E(b -y),V - M Vl 
and ages for these two clusters, although they may be physically associated (Vazquez et al. 1997). 

Hogg 22: The only saturated star in our images of Hogg 22 is HD 150958, an O type star and a 
speckle binary with angular separation of C'3 (McAlister et al. 1990). We find one possible Be star, 
No. 9, in Hogg 22. 

Hogg 22 is only 6' away from the cluster NGC 6204. While the two appear almost as the same group 
on the sky, Hogg 22 is the more distant yet less evolved cluster (Forbes & Short 1996). Because the 
two clusters share the same FOV in our images, we used the finder chart of Forbes & Short (1996) to 
consider the photometry of Hogg 22 members separately. However, the known members of Hogg 22 
do not include the full B-type MS, and as a result the cluster parameters have a large error. A more 
detailed analysis will be required to identify more cluster members. 

IC 2395: There are three saturated stars in our photometry of IC 2395: HD 74234, a B2 IV star; HD 
74455, a B1.5 Vn star; and HD 74531, a B2 V star in a double system. We identify only one possible 
Be star candidate in this cluster, No. 92. 

Our images of IC 2395 also contain members of the cluster BH 47, although both Claria et al. (2003) and 
Jorgcnsen & Westerlund (1988) suggest that IC 2395 and BH 47 are really a single open cluster within 
the Vela OB1C association. Therefore we consider the entire field without discriminating between the 
two clusters. 

IC 2581: The very bright, saturated star that dominates images of this cluster is V399 Car, a 
Cepheid variable. Two other bright stars, HD 90706 and HD 90707, were defocused slightly during our 
observations to avoid saturation. HD 90706 is a B1/B2 lab variable star (SIMBAD) and HD 90707 is 
a Bl III eclipsing binary (Lloyd Evans 1969). 

Our photometry led to the detection of two definite Be stars in IC 2581. The first, No. 158, is a B2 
Ve star, and the other is No. 243, a Bl Ve star. These two Be stars were also identified as Be stars 
by Lloyd Evans (1969, 1980). Lloyd Evans (1980) identifies No. 148 as a third Be star in IC 2581; 
however, it does not exhibit Ha emission in our results. We identify an additional 8 possible Be stars 
in IC 2581: Nos. 10, 123, 151, 157, 219, 267, 270, and 323. 

IC 2944: There are three saturated stars in our images of IC 2944. The first is HD 101131, an 06.5 
V((f)) + 08.5 V spectroscopic binary (Gies et al. 2002). In addition, HD 101205 is an 08 star and 
eclipsing binary of (3 Lyr type, and HD 101298 is an 06 V star. We identify No. 58 as a definite Be 
star in IC 2944, and Nos. 21 and 51 are possible Be stars. 

IC 2944 has a very distinctive appearance in our images; in addition to a slight amount of Ha nebulosity, 
several dark, sharp-edged Thackeray globules (Reipurth et al. 1997) are visible in the Ha images. 
Although these provide clear signs of active star formation in the field of IC 2944, there is disagreement 
within the literature about whether this is a significant physical cluster. While Walborn (1987) argues 
that the O stars of IC 2944 contribute to the ionization of an extensive HII region nearby, Perry & 
Landolt (1986) point out that both IC 2944 and IC 2948 are associated with this HII region in the 
Ccntaurus OB2 association, so the hot stars of IC 2944 are not required to ionize the region. They 
argue instead that IC 2944 is a chance superposition of early type field stars along the Carina arm 
of the Galaxy. Although the true nature of IC 2944 is not resolved, Figure 17 appears to represent 
a physically bound group, so we do include the cluster with our final statistics on Be stars in open 
clusters. 
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NGC 2343: No stars are saturated in our images of NGC 2343, and we do not identify any Be stars in 
this cluster. Although Claria (1972) identified our Nos. 19, 21, and 159 as probable foreground stars, 
they appear to lie on the cluster's MS and we classify these stars as B-type stars within the cluster. 
Also, he found that our No. 160 is probably not a cluster member, although we find it close to the MS 
and barely too red to classify as another B star. Finally, Claria (1972) identified No. 149 as a probable 
non-member of the cluster, but this star may be the only member of this sparse cluster to appear in the 
MS turnoff region, and confirmation of its membership in NGC 2343 could provide a better estimate 
of the cluster age. 

NGC 2362: Two stars in this cluster were saturated in our images, HD 57061 and HD 57192. HD 
57061 is a quadruple star system containing a visual double, two O-type stars with a separation of 
0'.'151. The brighter of the two O stars is a 154.9 day period spectroscopic binary that also contains 
an eclipsing binary with a 1.3 day period (van Leeuwen & van Gendcrcn 1997). The other saturated 
star in our images of NGC 2362 is HD 57192, a B2 V eclipsing binary (Kazarovets et al. 1999). We do 
not identify any definite Be stars in this cluster, and there is only one star, No. 19, that we classify as 
a possible detection. 

NGC 2367: There are no definite Be stars in this cluster, and only three possible candidates: Nos. 
55, 72, and 138. The two stars HD 57370 and BD -21°1881 are saturated in our images of NGC 2367. 
HD 57370 is a B2/B3 III star, and BD -21°1881 is a Bl V star. Although Ahumada & Lapasset 
(1995) classify HD 57370 as a blue straggler, we find several bright B-type stars in the cluster, and in 
fact NGC 2367 is one of the youngest clusters in this study. 

NGC 2383: Many of the stars in the nearby cluster NGC 2384 are in the same field of view as NGC 
2383, so we removed the 13 known members of the neighboring cluster from this set of photometry. 
The only saturated star in our images of NGC 2383 is the star HD 58509, but because it is a member 
of NGC 2384, its photometry is included with that cluster's data set. Two of the brightest stars in the 
field of NGC 2383 are not cluster members according to Subramaniam & Sagar (1999). They classify 
the first nonmember, No. 269, as an A3 I star. The other nonmembcr is No. 329, an Ml I star. Two 
definite Be stars are found in NGC 2383, Nos. 11 and 341. 

NGC 2384: Three known members of the nearby cluster NGC 2383 fell into the field of view of this 
cluster, so they were removed from the photometry set. None of the stars in NGC 2384 are saturated, 
and we found two possible Be stars in this cluster, Nos. 176 and 194. 

Although the two clusters NGC 2383 and NGC 2384 are separated by only about 5', they are not 
associated. Vogt & Moffat (1972) and Subramaniam & Sagar (1999) have shown that NGC 2384 is 
both younger and at a farther distance than its neighbor, indicating a chance superposition in the plane 
of the Galaxy. According to Vogt & Moffat (1972), the two clusters have nearly the same reddening 
because NGC 2383 lies on the outer edge of the local spiral arm, so there is little dust between it and 
the more distant NGC 2384. 

NGC 2414: The star HD 60308, a B2 lab star, is saturated in our images of NGC 2414. There are 
five definite Be stars in this cluster (Nos. 12, 65, 73, 87, and 97). In addition, we found 13 possible 
Be star candidates (Nos. 2, 3, 4, 14, 17, 19, 34, 37, 40, 106, 114, 119, and 122). However, there are a 
large number of apparent foreground stars contaminating the MS of this cluster, so we doubt that the 
Be candidates are true Ha emitters. 

In fact, we doubt that this is a true, physically associated cluster. Although the group appears 
somewhat prominent against the field, we see two possible main-sequences in the color-magnitude 
diagram that we attempt to fit with different values for E(b — y) and Vb — My. Neither value of 
E(b — y) provides a good match between the observed color-color sequence and the theoretical color- 
color curve, which suggests that the group either has a highly variable reddening or is not associated. 
We see no evidence for nebulosity or dust gradients in our Ha images of the cluster. 
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NGC 2421: There are four stars in NGC 2421 with definite Ha emission and b — y < 0.350, the cutoff 
for OB stars based on the cluster's reddening. These definte Be stars are Nos. 311, 356, 378, and 396. 
Slettebak (1985) observed H7 emission from No. 356 and classified its spectral type as B3— B5 (V:). 
He also observed weak H7 central emission in No. 311, which he also classified as a B3— B5 (V:) star. 
Also, Stephenson & Sanduleak (1977) identified No. 396 as a B9 star with weak emission lines. No. 
378 is a newly identified Be star. 

We also find six possible Be stars in NGC 2421: Nos. 83, 197, 256, 341, 483, and 530. Stephenson 
& Sanduleak (1977) also observed weak emission in No. 341. In addition, Stephenson & Sanduleak 
(1977) and Slettebak (1985) saw weak emission in No. 264, although this star is both too red and too 
faint in Ha to be a Be star in our results. 

No stars in our images of NGC 2421 were saturated. 

NGC 2439: In our images of NGC 2439, there were two stars that had to be defocused to avoid 
saturation. The brightest is R Pup, an F8/G0 la star. Eggen (1983) describes R Pup as a "pseudo- 
cepheid" , a non-periodic variable supergiant whose low dispersion spectrum resembles a Cepheid. The 
second brightest star in the cluster is V384 Pup, an M3 IaO-Ia pulsating variable star. 

Our photometry of NGC 2439 revealed six definite Be stars, Nos. 22, 24, 58, 101, 193, and 196. No. 
24 was classified as a Bl V: star with weak central emission in the H7 line by Slettebak (1985). He 
also observed double peaked Ha emission and sharp central absorption cores in the H7 and B.S lines 
of No. 58. Slettebak (1985) also noted that the star No. 101 had been previously identified as a Be 
star, presumably in an unpublished catalog by Sanduleak & Robertson (1977). According to Slettebak, 
Sanduleak & Robertson list a total of five Be stars in NGC 2439, although the remaining two arc not 
identified by Slettebak. In addition to the six definite Be stars found from our photometry, we detect 
seven possible cases: Nos. 13, 16, 43, 91, 116, 147, and 199. 

NGC 2439 may not be a true open cluster (Eggen 1983; Kaltcheva, Gredel, & Fabricius 2001). 
Kaltchcva ct al. (2001) found many localized clumps of OB stars within the cluster as well as highly 
variable reddening. In addition, they found a large scatter in distances, proper motions, and radial ve- 
locities among their sample of early type stars in the field. While NGC 2439 may simply be a collection 
of field stars revealed by an absorption hole in the galactic plane, its color- magnitude and color-color 
diagrams in Figure 25 appear to represent a true physical cluster, so we do include the cluster in our 
analysis in §3 and §4. 

NGC 2483: There were no stars that were saturated in our shortest exposures of this cluster, so 
there is no chance of any Be stars being overlooked. One of the definite Be stars in this cluster, No. 
186, was classified as an AO I emission line star by Stephenson & Sanduleak (1977), but they noted an 
uncertainty in the spectral type and a Be classification is likely. In addition, we identify No. 87 as a 
Be star, and there are ten possible Be stars in this cluster: Nos. 7, 27, 40, 42, 53, 90, 97, 179, 290, and 
292. 

NGC 2483 is not a true, physically associated cluster, and we ignore it in our analysis. FitzGcrald 
& Moffat (1975) noted that the region has only a slight increase in the number density of stars over 
the background field, and they identified three distinct groups of stars, each at different distances, 
that are superimposed on the field. Havlen (1976) confirmed these results and recommended that 
NGC 2483 no longer be classified as a galactic cluster. Havlen also showed that the four bluest stars 
in our photometry are a group of foreground B6 — B8 stars that are somewhat less reddened (mean 
E(B — V) = 0.11) than the population of 9 09 — B5 stars that we used to fit the "cluster's" isochrone. 

NGC 2489: There are no Be stars in this cluster, and no stars are saturated. 

NGC 2571: In our images of NGC 2571, no stars were saturated. The single Be star in this cluster 
is No. 76, and it is classified as a B4 star in SIMBAD. 
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NGC 2659: None of the stars in NGC 2659 were saturated in our images. One definite Be star 
emerged from our color-color plot of this cluster, No. 233. In addition to this Be star, there are 4 
possible Be stars (Nos. 3, 42, 89, and 222). 

NGC 3293: NGC 3293 is a young cluster with ongoing star formation around the cluster's periphery 
(Baume et al. 2003), although its nuclear age is used in this study. Baume et al. found many pre-MS 
members in this cluster, as well as several massive stars that have already begun to move off of the 
ZAMS. Both they and Turner et al. (1980) observed variable reddening in this cluster, and our own 
Ha image of NGC 3293 reveals thin, wispy clouds of gas across the region. The differential reddening 
is a likely cause for error in the isochrone fit, and the background Ha emission as well as emission from 
pre-MS stars is also likely to contaminate the number of Be stars detected in this survey. 

To further amplify the problem of detecting Be stars in NGC 3293, a large number of stars are saturated 
in our images of this bright cluster. They are: HD 91943, a B0.5/B1 II/III variable star; HD 91969, a 
BO lb star; HD 92007, a Bl III variable star of A Eri type and a probable Be star (Balona 1994); HD 
92044, a Bl II star; CD -57°3348, a B0.5 V star; CPD -57°3506B, a Bl III star that is a suspected 
eclipsing binary and a variable star of (3 Cep type (Balona 1994); CPD — 57°3521, a Bl III star; 
CPD -57°3523, a Bl III star; CPD -57°3526B, a B0. 5 III star; and V361 Car, an M0 star. Clearly, 
the Be star population is probably not well sampled in this cluster due to the large number of B-type 
giants that are omitted from the photometry, including at least one probable Be star. 

We find three definite Be stars in NGC 3293: Nos. 6, 42, and 46. In addition, there are eight possible 
Be stars with somewhat weaker emission: Nos. 4, 8, 27, 51, 54, 84, 99, and 116. No. 99 was noted to 
be a Be star by Shobbrook (1980), and No. 51 was also identified as a Be star by Feast (1958). 

NGC 3766: The cluster NGC 3766 was the test case for our photometric search for Be stars, and it 
was discussed in detail in Paper 1. The five definite Be stars that we identified with our photometry 
are (in our numbering scheme) Nos. 47, 127, 154, 198, and 200. Our spectra confirmed these five and 
also identified three other Be stars with weak Ha emission: Nos. 31, 83, and 92. From our photometry, 
we classified Nos. 61, 73, 83, 92, 94, 126, 130, and 139 as possible Be detections. Our spectroscopy 
proves that at least some of these possible Be stars are indeed weak emission line stars. 

The five saturated stars in our images of NGC 3766 are V910 Con, HD 100840, HD 100943, HD 306794, 
and HD 306799. One of these saturated stars, HD 100943, was identified as a Be star by Shobbrook 
(1987) and confirmed by our spectra. However, for consistency in the detection technique, we only 
include Be stars detected with our photometry in this study. 

NGC 4103: No stars in our images of NGC 4103 were saturated. One star, No. 231, clearly exhibits 
Ha emission in the color-color diagram of this cluster. It was classified as a B2 IVe star by Wessclink 
(1969), and we classify it as a definite Be star detection. Two other possible Be stars (Nos. 13 and 
209) were also found in this cluster. 

NGC 4755: In our images of NGC 4755, six bright stars are saturated. They are: HD 111904, a B9 
la variable star; HD 111934, a B2 lb variable star; HD 111973, a B5 la star; HD 111990, a B1/B2 
lb star in an optical double system; SAO 252073, an M2 lab variable star of irregular type; and SAO 
252075, a B0. 5 Vn ellipsoidal variable star. 

Of the remaining unsaturated stars, five are definitely Be stars. They are: Nos. 43, 129 (B1.5 Vpne; 
Schild 1970), 171 (B2 IVnc; Knoechel 1980), 180, and 187. Our photometry detected 7 additional Be 
star candidates (Nos. 20, 53, 131, 168, 189, 222, and 223). 

NGC 5281: One star, HD 119699, is saturated in our images of this cluster, and it is classified as an 
Al II star. The only definite Be star detected in our study is No. 74, and five additional Be candidates 
were found: Nos. 11, 13, 29, 66, and 79. The bright star No. 80 also exhibits Ha emission that 
distinguishes it in the color-color diagram, and it was previously classified as a Be star by Humphreys 
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(1975) and as a blue straggler by Ahumada & Lapasset (1995). However, its (b — y) color places it 
slightly redder than the B-type stars in our study so we do not classify it as a Be star, nor does it 
appear to be a blue straggler in the cluster. 

NGC 5593: There are no saturated stars in our images of NGC 5593. The star No. 204 is a definite 
Be star in this cluster, and No. 180 is a possible Be star. 

NGC 6178: The only saturated star in our images of NGC 6178 is HD 149277, a B2 IV/V star, and 
we do not identify any Be stars in this cluster. 

NGC 6193: This young cluster is rich in spectroscopic binaries, and Arnal et al. (1988) found a binary 
fraction as high as 72%. Two of the saturated stars in NGC 6193 are HD 150135 and HD 150136, an 
06.5 V + 05 visual binary and the only two O stars in the Ara OBI association (Hcrbst & Havlen 
1977). In addition, the B1/B2 lb supergiant HD 150041 is saturated in our images. NGC 6193 is an 
extremely young cluster, about the same age as the Orion nebula cluster (Herbst & Havlen 1977), with 
pre-MS B stars moving onto the ZAMS at magnitudes V > 12.4 (Moffat & Vogt 1973). Therefore the 
possible Be star candidate in NGC 6193, No. 38, may be a pre-MS Herbig Be star. 

NGC 6200: Although the cluster NGC 6200 barely stands out against the background field, it contains 
a large number of B-type stars and Be star candidates. We identify Nos. 2, 38, 112, 213, 239, 247, 259, 
and 305 as possible Be stars and Nos. 115 and 153 as definite Ha emission stars. NGC 6200 appears 
close to NGC 6204 and Hogg 22 on the sky, although this is a chance superposition since the cluster is 
younger and more distant than its apparent neighbors. 

NGC 6204: NGC 6204 is only 6' away from Hogg 22 on the sky (see the notes for that cluster). 
Members of Hogg 22 were removed from our results for NGC 6204, and no members of this cluster are 
saturated in our images. While we identify no definite Be stars in NGC 6204, there are two possible 
Be stars: Nos. 46 and 47. No. 4 was also identified as an Ha emitter by Vega et al. (1980), but it is 
barely too red to be a B-type star according to our results. 

NGC 6231: NGC 6231 is a very bright cluster, and a large number of massive stars are saturated in 
our images. They are: HD 152233, an 06 III variable star; HD 152234, a B0.5 la star; HD 152248, 
an 07 I + 07 I spectroscopic binary (Penny, Gies, & Bagnuolo 1999); HD 152270 (WR 79), a WC+ 
Wolf-Rayet star; HD 326331, a B star; HD 152314, an 0+ variable star; HD 152219, an 09 IV 
spectroscopic binary; SAO 227385, an 09 III spectroscopic binary; HD 152249, an 09 lb supergiant; 
and HD 152218, an 09 V spectroscopic binary. The cluster has an extremely high binary frequency; 
Garcia & Mermilliod (2001) found that 82% of the members earlier than B1.5 V are binary. It also 
contains a large number of pulsating variable stars, and Arentoft et al. (2001) identify at least 14 stars 
with (3 Cephei, 6 Scuti, 7 Doradus, and other types of pulsations. We identify two definite Be stars, 
Nos. 67 and 79, and three possible Be stars, Nos. 2, 35, and 221, but none of these are identified as 
pulsators by Arentoft et al. (2001). Sung, Bessell, & Lee (1998) found 19 pre-MS stars or candidates 
in NGC 6231, including our Be star No. 79. 

NGC 6249: No stars are saturated in our images on NGC 6249, and we identify one definite Be star, 
No. 207. 

NGC 6250: There are two saturated stars in our images of NGC 6250. They are HD 152917, an A6 
III/IV star, and HD 152853, a B2 III star. The cluster is sparse, with only seven B-type stars, and we 
find no Be stars in NGC 6250. 

NGC 6268: NGC 6268 is an evolved cluster with very few B stars remaining. We identify one of 
them, No. 51, as a definite Be star. The only saturated star in our images of this cluter is HD 322550, 
an FO star and likely foreground to the cluster. 
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NGC 6322: Three stars in NGC 6322 are saturated in our images, HD 156292, a Bl Ib/II supergiant, 
HD 156234, a BO III star, and HD 156189, an A3 III/IV star. The definite Be star in this cluster is No. 
12, and Nos. 7, 26, 50, 51, 76, 79, 84, and 85 are possible candidates. Most of these are likely Herbig 
Be stars as pre-MS contraction is still occuring in NGC 6322 among stars with V > 12.5 (Moffat & 
Vogt 1975c). 

NGC 6425: The only saturated star in our images of NGC 6425 is HD 161508, a G6/G8 V star and 
likely a foreground star. We find no Be stars in this cluster. 

NGC 6530: NGC 6530 is part of a complex star-forming region, the Lagoon Nebula, and many of the 
bright B stars in the cluster are pre-MS stars moving onto the ZAMS (Sung, Chun, & Bessell 2000). 
The B0 IVpne star HD 164906 (Hiltner, Morgan, & Neff 1965) is saturated in our images of NGC 6530. 
We find one definite Be star, No. 50, and two possible Be stars, Nos. 20 and 31, in the cluster. No. 
50 was classified as a B2 Ve star by Hiltner et al. (1965) and it was identified as a Herbig Be star by 
Boesono, The, & Tjin A Djie (1987). 

NGC 6531: The two stars HD 164863 and HD 164844 are saturated in our images of NGC 6531. 
They are 0+ and B1/B2 III stars, respectively. Although Schild & Romanishin (1976) found no Be 
stars in this cluster, we see Ha emission from No. 12. Even though we classify this as a definite Be 
star because its y — Ha color is more than 5 standard deviations above the other B star colors, the 
photometry errors for this faint star are comparable to the strength of the emission so it may be a false 
Be star detection. 

NGC 6604: NGC 6604 appears faint without a well-defined MS in our images due to its location at 
the core of the HII region S54 (Barbon et al. 2000). The saturated star in our images is HD 167971, 
an 05—8 V + 05—8 V + 08 lb triple system (Lcithcrcr et al. 1987) and eclipsing binary of [3 Lyr 
type. The cluster may have some faint pre-MS members with V > 16 (Barbon et al. 2000). There are 
seven possible Be stars in NGC 6604, Nos. 2, 3, 9, 11, 21, 66, and 73, although they all lie at the faint 
end of the MS and the errors in their y — Ha colors are relatively high. We believe that the apparent 
emission is more likely due to photometric scatter. 

NGC 6613: NGC 6613 = M18 is a sparse cluster in Sagittarius. No stars are saturated in our images 
of this cluster, and we find one definite Be star, No. 1, also discovered as an Ha emission star by 
MacConnell (1981). 

NGC 6664: No stars are saturated in our images of NGC 6664, and we find several Be star candidates. 
Nos. 177 and 221 have definite Ha emission, and Nos. 137 and 182 are possibly Be stars as well. Star 
No. 176 also appears to have Ha emission, but it is not classified as a B-type star in our results. 

Ruprecht 79: There is one star that is barely saturated in our b images of Ruprecht 79, but not 
much is known about it. The SIMBAD database identifies it as GSC 08585-01823 with B = 10.9 and 
V = 10.28. Its color, B — V = 0.6, suggests that it is either a massive OB star in the cluster or an 
unreddened, foreground star. The latter is more likely, given the presence of many stars to the left of 
the MS which must be foreground stars. Walker (1987) points out that unreddened G and K stars fall 
into the same region as the reddened B and A stars in his color-color diagram of this cluster. 

In a color-magnitude diagram of Ruprecht 79, the MS is not clearly defined due to the overlap of 
many foreground G and K stars (Walker 1987). They also contaminate the color-color plot, making 
it difficult to fit a theoretical reddening curve to our data. Because of the contamination, there are 
a large number of stars that lie more than 2a above the theoretical reddening curve in the Be star 
region of the color-color plot, but they are not likely to be Be stars since there is no clear gap in the 
stellar population indicating a group of Ha emitters. Therefore the Be star detections in this cluster 
are likely incorrect, and we disregard these detections in the analysis below. 
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The cluster Ruprecht 79 has been studied a handful of times in the literature (Moffat & Vogt 1975a; 
Harris & van den Bergh 1976; Walker 1987) because the Cepheid variable CS Velorum is a member of 
the cluster. However, Harris & van den Bergh (1976) note that the cluster is "rather loose", and they 
point out that no large-scale photographs of the region were available at the time of their observing 
program. That is unfortunate, because a 60' x 60' image of Ruprecht 79 from the STScI Digitized 
Sky Survey reveals that Ruprecht 79 is not at all prominent against the field. In fact, the eastern half 
of the image is more crowded with stars than the west, suggesting a gradient in the dust absorption 
across the field. From our eyeball estimates, there are many other regions in the field with comparable 
stellar density to this cluster. 

In addition, the galactic coordinates of Ruprecht 79 and its distance (see Table 4) place it within the 
Sagittarius-Carina spiral arm of the Milky Way (Binney & Mcrrificld 1998). The cluster is viewed 
through the tangential edge of the arm (Dame et al. 1987; Binney & Merrifield 1998), so a gradient in 
absorption across the field of view is likely Therefore we conclude that Ruprecht 79 simply corresponds 
to a hole in the dust in this region of the Galaxy, and it is not a true open cluster. 

Ruprecht 119: Although Moffat & Vogt (1973) classified Ruprecht 119 as a non-cluster, Piatti et al. 
(2000b) showed that it is in fact physically bound. We find a distinct MS for the cluster, although we 
confirm a slight gap in the MS around y — 13.5, and possibly at y — 16, both identified by Piatti et 
al. (2000b). No stars are saturated in our results for Ruprecht 119, and we find one definite Be star 
(No. 332) and nine possible Be stars (Nos. 267, 301, 412, 427, 487, 549, 574, 674, and 706). 

Ruprecht 127: There are no saturated stars in our photometry of Ruprecht 127. We find no clear 
MS visible from our data, and any true MS is somewhat contaminated by field stars. We experimented 
with fitting the MS using the range of E(b — y), Vq — My, and ages available in the literature (see Table 
2) as well as our own values for a less reddened, closer group of stars. Each isochrone fit illustrated that 
multiple, sparse stellar populations are in the FOV, and there appear to be two stellar populations, 
with distinctly different reddening, in the color-color sequence as well. The Be star candidates found 
using the parameters in Table 2 are likely foreground stars rather than true emission stars. We find no 
Be stars when we fit the more nearby population. We doubt that this is a physical cluster, and instead 
we believe it is a chance physical alignment of several bright stars. 

Ruprecht 140: Only the star GSC 07397-00784 is saturated in our images of Ruprecht 140. We find 
no clear MS in this cluster, and we agree with Piatti & Claria (2001) who concluded that Ruprecht 
140 is not a physical cluster but rather a density fluctuation over the background field. We find one 
possible Be star in this FOV (No. 468). 

Stock 13: We find one definite Be star, No. 831 , and one possible Be star, No. 271, in our results from 
Stock 13. No stars are saturated in our data. There may be some pre-MS stars in this cluster, as Moffat 
& Vogt (1975b) identified both Nos. 376 and 452 as possible pre-MS stars. We see a distinct cluster 
MS, although the color-color diagram possibly indicates a small population of nearly unreddened, late 
B- and early A- type foreground stars that are also present. Steppe (1977) also found a group of 20 
foreground red giants in the field of Stock 13. The presence of foreground stars is not surprising 
considering the relatively large distance to this cluster (see Table 4). 

Stock 14: In our images of Stock 14, the F8 la supergiant V810 Cen is saturated, although Kienzle 
et al. (1998) found that this star is not a member of the cluster. We find one definite Be star (No. 
164) and a possible Be star (No. 105) in Stock 14, and both of these stars are also identified as Be or 
emission line stars in SIMBAD. 

Trumpler 7: No definite Be stars are found in this cluster, but 5 possible Be stars have Ha emission 
more than 2a above the color-color fit. These Be star candidates are Nos. 13, 55, 57, 59, and 82. Only 
the Al IV star HD 59163 is saturated in our images of Trumpler 7. 
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Trumpler 18: No stars are saturated in our results from Trumpler 18, and we identify five possible 
Be stars (Nos. 10, 29, 55, 137, and 219). Like Vazquez & Feinstein (1990), we see a clear MS for 
this cluster. However, Vazquez & Feinstein (1990) questioned the physical association of Trumpler 
18 because they found a peculiar luminosity function for the cluster. We compute the mass function 
for Trumpler 18 based on the method described in §3, and we find a normal power law relation with 
T = —1.5 ± 0.16, depending on the number of logM bins used. With this value of T we should expect 
75 ± 4 B-type stars in the cluster, and we detect 71 (see Table 2). We conclude that Trumpler 18 is a 
true open cluster. 

Trumpler 20: There are no saturated stars in our images of Trumpler 20. We find very few B-type 
stars in this cluster, and only one possible Be star (No. 1039). 

Trumpler 27: We find two saturated stars in our images of Trumpler 27: HD 318014, a B8 lab 
supergiant, and CD — 33°12241, an M0 la binary. The cluster contains one star with very strong 
emission, No. 15 (WR 95). In this case, the emission more likely originates from the He II A 6527 line 
than Ha in this Wolf-Rayet star. We also find two possible Be stars, Nos. 49 and 52. We do not see a 
well defined MS for Trumpler 27, although we believe this is due to the cluster's significant differential 
reddening. The & Stokes (1970) found that the cluster lies behind the dust lane B217 that partially 
encircles the cluster, and Feinstein et al. (2000) used polarization measurements to determine that two 
separate dust components lie in the field of Trumpler 27. We believe that our estimates of the cluster 
parameters are unreliable, and we disregard this cluster in our statistical analysis of Be stars (§3). 

Trumpler 28: The B3 star HD 317971 is saturated in our images of Trumpler 28, and we find no Be 
stars in this cluster. 

The cluster Trumpler 28 is only 14' from NGC 6383, and the two may form a binary open cluster system 
(Subramaniam et al. 1995). However, the large angular size of NGC 6383 (20'; WEBDA) excluded 
it from our target list, and we do not present photometry of that cluster in this study. According 
to Subramaniam et al. (1995), NGC 6383 is significantly younger than Trumpler 28, and the three 
apparent blue stragglers near the MS of Trumpler 28 (Nos. 139, 161, and 193) are probably members 
of NGC 6383 instead. 

Trumpler 34: No stars are saturated in our images of Trumpler 34. We find one definite Be star, 
No. 355, and one possible Be star, No. 297, in this cluster. No. 355 is very blue for a late B-type star, 
but its photometric errors are small enough that we trust its position in Figure 62. However, it may 
possibly be a foreground star. 

vdB-Hagen 217: There are no saturated stars in our images of vdB-Hagen 217. We find no Be stars 
in this cluster. 
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Fig. 1 . — The distribution of the true open clusters 
is shown relative to the Sun (©) and the Galac- 
tic center (x). The positions of the spiral arms 
(gray) are from Cordes & Lazio (2003), and the 
orientation of the central bar (dotted line) is from 
Bissantz, Englmaier, & Gerhard (2003). 




Fig. 2. — The percentage of Be stars (diamonds) 
is shown MS il function of their distance from the 
Galactic center, d cen t- The upper limit for each 
cluster, the percentage of definite and possible Be 
stars, is also shown. 
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Fig. 3. — The percentage of Be stars is shown as a 
function of the cluster density in the same format 
as Figure 2. 
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Fig. 5. — The average definite Be frequency (solid 
lines) and total possible Be frequency (dashed 
lines) are plotted for four age bins correspond- 
ing to log age < 7.0, 7.0 < log age < 7.4, 7.4 < 
log age < 8.0, and log age > 8.0. The errors in Be 
frequency represent the Poisson statistics for our 
sample, and the age error bars represent the max- 
imum and minimum ages observed in each bin. 
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Fig. 4. — The percentage of definite Be stars is 
shown as a function of cluster age in the same 
format as Figure 2. 
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Fig. 6. — The definite Be star detections are shown 
on an absolute color-magnitude diagram. Three 
isochrones with ages 10 5 , 10 7 , and 10 8 yr are 
shown for comparison. The B stars are divided 
into four regions: evolved, early, mid, and late 
spectral types. The contours represent the dis- 
tribution of normal B-type stars detected in our 
sample (29%, 43%, 57%, and 76%, respectively). 
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Fig. 7. — A histogram plot of the total Be and B 
star populations of all clusters. The bins are di- 
vided into evolved, early, mid, and late B spectral 
types. 
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Fig. 8. — My of the definite Be stars is plotted as 
a function of cluster age (diamonds). The abso- 
lute magnitude of a B star at the TAMS is also 
shown (dashed curve; Lejeune & Schaerer 2001). 
For a primary at the TAMS and system mass ratio 
q = 0.3, My of the secondary is also shown (dotted 
curve). 




Fig. 9. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Basel 1. The 
isochrone fit (solid line in top), color-color fit 
(solid line in bottom), reddening value (dotted 
lines), yBmin and ysmax (dashed lines in top), and 
parabolic fit of unrcddened stars in the atlas of Ja- 
coby et al. (1984) (dashed line in bottom) are also 
shown. Be stars (large) are distinguished from all 
other stars (small) in both diagrams (although in 
this cluster no Be stars are present). Potential 
foreground stars or white dwarfs are also marked 
( x ) , although none are present in our data for this 
cluster. 
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Fig. 10. — Color-magnitude {top) and color-color 
( bottom) diagrams of the cluster Bochum 13 in the 
same format as Figure 9. 



Fig. 12. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster Haffncr 16 in the 
same format as Figure 9. 




Fig. 11. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Collinder 272 in 
the same format as Figure 9. 



Fig. 13. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Hogg 16 in the 
same format as Figure 9. 
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Fig. 14. — Color-magnitude {top) and color-color 
(bottom) diagrams ol the cluster Hogg 22 in the 
same format as Figure 9. 



Fig. 16. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster IC 2581 in the 
same format as Figure 9. 




Fig. 15. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster IC 2395 in the 
same format as Figure 9. 



Fig. 17. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster IC 2944 in the 
same format as Figure 9. 
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Fig. 18. — Color-magnitude {top) and color-color 
(bottom) diagrams ol the cluster NGC 2343 in the 
same format as Figure 9. 



Fig. 20. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 2367 in the 
same format as Figure 9. 




Fig. 19. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2362 in the 
same format as Figure 9. 



Fig. 21. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2383 in the 
same format as Figure 9. 
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Fig. 22. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 2384 in the 
same format as Figure 9. 



Fig. 24. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 2421 in the 
same format as Figure 9. 




Fig. 23. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2414 in the 
same format as Figure 9. 



Fig. 25. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2439 in the 
same format as Figure 9. 
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Fig. 26. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2483 in the 
same format as Figure 9. 



Fig. 28. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2571 in the 
same format as Figure 9. 




Fig. 27. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2489 in the 
same format as Figure 9. 



Fig. 29. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 2659 in the 
same format as Figure 9. 
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Fig. 30. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 3293 in the 
same format as Figure 9. 



Fig. 32. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 4103 in the 
same format as Figure 9. 




Fig. 31. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 3766 in the 
same format as Figure 9. 



Fig. 33. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 4755 in the 
same format as Figure 9. 
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Fig. 34. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 5281 in the 
same format as Figure 9. 



Fig. 36. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6178 in the 
same format as Figure 9. 




Fig. 35. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 5593 in the 
same format as Figure 9. 



Fig. 37. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6193 in the 
same format as Figure 9. 
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Fig. 38. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 6200 in the 
same format as Figure 9. 



Fig. 40. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 6231 in the 
same format as Figure 9. 




Fig. 39. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6204 in the 
same format as Figure 9. 



Fig. 41. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6249 in the 
same format as Figure 9. 



31 




Fig. 42. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 6250 in the 
same format as Figure 9. 



Fig. 44. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster NGC 6322 in the 
same format as Figure 9. 




Fig. 43. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6268 in the 
same format as Figure 9. 



Fig. 45. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6425 in the 
same format as Figure 9. 
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Fig. 46. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6530 in the 
same format as Figure 9. 



Fig. 48. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6604 in the 
same format as Figure 9. 




Fig. 47. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6531 in the 
same format as Figure 9. 



Fig. 49. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6613 in the 
same format as Figure 9. 
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Fig. 50. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster NGC 6664 in the 
same format as Figure 9. 
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Fig. 51. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Ruprecht 79 in 
the same format as Figure 9. 
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Fig. 52. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Ruprecht 119 in 
the same format as Figure 9. 
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Fig. 53. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Ruprecht 127 in 
the same format as Figure 9. 
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Fig. 54. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Ruprecht 140 in 
the same format as Figure 9. 



Fig. 56. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Stock 14 in the 
same format as Figure 9. 




Fig. 55. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Stock 13 in the 
same format as Figure 9. 



Fig. 57. — Color-magnitude (top) and color-color 
( bottom) diagrams of the cluster Trumplcr 7 in the 
same format as Figure 9. 
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Fig. 58. — Color-magnitude {top) and color-color 
(bottom) diagrams of the cluster Trumpler 18 in 
the same format as Figure 9. 



Fig. 60. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Trumpler 27 in 
the same format as Figure 9. 




Fig. 59. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Trumpler 20 in 
the same format as Figure 9. 



Fig. 61. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Trumpler 28 in 
the same format as Figure 9. 
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Fig. 62. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster Trumpler 34 in 
the same format as Figure 9. 




Fig. 63. — Color-magnitude (top) and color-color 
(bottom) diagrams of the cluster vdB-Hagen 217 
in the same format as Figure 9. 



Table 1 

Photometry of Open Clusters 



Cluster 


No. 


RA (2000) 


Dec (2000) 


V 


Sy 


b-y 


S(b-y) 


y — Ha 


5(y - Ha) 


Code WEBDA Identifier 


Bascll 


1 


18 47 56.08 


-5 45 51.0 


16.445 


0.066 


0.988 


0.103 


0.917 


0.103 


O 


Basell 


2 


18 48 01.11 


-5 45 50.7 


15.981 


0.056 


0.744 


0.084 


0.770 


0.094 


o 


Bascll 


3 


18 48 23.49 


-5 45 50.9 


16.137 


0.060 


0.868 


0.092 


0.876 


0.098 


o 


Bascll 


4 


18 48 33.34 


-5 45 53.7 


16.570 


0.069 


0.770 


0.102 


0.898 


0.114 


o 


Bascll 


5 


18 48 06.99 


-5 45 55.7 


14.437 


0.050 


1.396 


0.079 


1.537 


0.064 


o 


Bascll 


6 


18 48 37.87 


-5 45 53.8 


15.479 


0.048 


0.657 


0.070 


0.693 


0.081 


o 


Basell 


7 


18 47 52.85 


-5 45 56.7 


16.315 


0.064 


0.868 


0.096 


1.037 


0.096 


o 


Bascll 


8 


18 48 01.30 


-5 45 57.1 


14.802 


0.047 


1.111 


0.074 


1.181 


0.066 


o 


Bascll 


9 


18 48 31.17 


-5 45 56.8 


15.899 


0.054 


0.812 


0.084 


0.683 


0.093 


o 


Basell 


10 


18 48 34.64 


-5 45 56.7 


13.519 


0.031 


0.429 


0.046 


0.395 


0.052 


o 



NOTE. — The complete version of this table is in the electronic edition of the Journal. For each star observed, we give the cluster name and our 
identifying number, the right ascension (RA), and declination (Dec) for the epoch 2000. We also provide the y magnitude, the b — y and y — Ha colors, 
and the error for each. The code is used to label definite Be stars (Be), possible Be stars (Be?), B-typc stars (B), probable foreground stars (F), and 
other stars in the field (O). Finally, the WEBDA number and other identifiers arc provided where available. 



Table 2 
Summary of Open Clusters 



Cluster 


RA (2000) 


Dec (2000) 


E(b 


ail 

y> 




1V1 \r 


log 


Age 


Mr, 


Mi-, 

(Possible) 


-< v B + Be 


Ixcicrcnccs 


Basel 1 


18 


48 


12 


-05 


51 


00 





33 


11 


00 


8 


45 








3 


6 


Bochum 13 a . . 


17 


17 


24 


-35 


33 


00 





64 


10 


16 


6 


82 


2 


l 


> 23 


1 


Collindcr 272 . 


13 


30 


26 


-61 


19 


00 





34 


11 


55 


7 


11 


3 


5 


91 


1,30 


Haffncr 16 .... 


07 


50 


20 


-25 


28 


00 





13 


12 


50 


7 


08 





1 


19 


1,7,14 


Hogg 16 


13 


29 


18 


-61 


12 


00 





31 


11 


00 


7 


05 





1 


35 


1 


Hogg 22 


16 


46 


37 


-47 


05 


00 





47 


11 


75 


6 


70 





1 


> 22 


8,14 


IC 2395 


08 


42 


37 


-48 


06 


48 





05 


9 


24 


7 


22 





1 


14 


1 


IC 2581 


10 


27 


24 


-57 


38 


00 





31 


12 


00 


7 


14 


2 


8 


117 


1,13,14 


IC 2944 b 


11 


38 


20 


-63 


22 


22 





24 


11 


27 


6 


82 


1 


2 


50 


1 


NGC 2343 


07 


08 


18 


-10 


39 


00 


o 


11 


10 


12 


8 


00 


o 





11 


1,11,14 


NGC 2362 


07 


18 


48 


-24 


57 


00 


{) 


07 


10 


71 


6 


91 


o 


1 


41 


1 


NGC 2367 


07 


20 


06 


-21 


53 


00 





20 


11 


84 


6 


74 


o 


3 


16 


1,7,31 


NGC 2383 


07 


24 


41 


-20 


56 


42 


{) 


20 


12 


62 


7 


40 


2 





20 


14,26,31 


NGC 2384 


07 


25 


12 


-21 


01 


24 


() 


22 


11 


63 


6 


90 





2 


33 


1,31 


MPP 941/l b 


07 


33 


XI 


— 15 


27 


12 





41 


1 1 


99 


g 


80 


5 


13 


86 


14 31 


wpp 9491 


07 


36 


18 


— 20 


37 


00 




35 


1 1 


69 


7 


40 




g 


65 


1 14 16 




u / 


i n 
411 


4o 


— 31 


39 


00 





30 


12 


18 


7 


25 


6 


7 


111 


1,14 


NGC 2483 b . . . 


07 


55 


54 


— 27 


56 


00 





30 


12 


31 


7 


00 


2 


10 


49 


14 


NGC 2489 .... 


07 


56 


12 


—30 


04 


00 





30 


10 


81 


8 


38 








7 


14,19 


NGC 2571 


08 


18 


04 


— 29 


A A 

44 


00 





10 


10 


64 


7 


49 


1 





24 


1 


NGC 2659 


08 


42 


36 


—44 


57 


00 





38 


10 


80 


6 


89 


1 


4 


32 


1,27 


NGC 3293 .... 


10 


35 


49 


— 58 


13 


00 





22 


12 


20 


6 


90 


3 


8 


> 91 


4 


IN GO O/DO .... 


11 


36 


13 


-61 


36 


55 





15 


1 1 


40 


7 


16 


5 


8 


> 146 


1,25 


L\kj\^i 41UO .... 


12 


06 


43 


-61 


15 


21 


II 


22 


X i 


i U 


"7 

i 


•50 


1 


9 


Q9 




1 N vj V i yu .... 


12 


53 


42 


— 60 


22 


00 


(} 


28 


11 


80 


g 


85 


5 


7 


> 178 


3 


NGC 5281 


13 


46 


30 


-62 


54 


54 





19 


11 


00 


7 


71 


1 


5 


40 


14,15 


NGC 5593 


14 


25 


54 


-54 


49 


00 





26 


10 


62 


7 


00 


1 


1 


15 


34 


NGC 6178 


16 


35 


42 


-45 


38 


00 





18 


9 


80 


7 


60 





{) 


6 


15,21 


NGC 6193° . . . 


16 


41 


18 


-48 


46 


00 





33 


10 


31 


6 


48 





1 


31 


1,10,15 


NGC 6200 


16 


44 


12 


-47 


29 


00 





43 


11 


56 


6 


93 


2 


8 


112 


1 


NGC 6204 


16 


46 


08 


-47 


00 


44 


o 


35 


10 


40 


8 


10 





2 


34 


5,8 


NGC 6231° 


16 


54 


09 


-41 


49 


36 


o 


33 


10 


47 


6 


90 


2 


3 


129 


1,20 


NGC 6249 


16 


57 


36 


-44 


47 


00 


o 


34 


10 


06 


7 


40 


1 





13 


14,15 


NGC 6250 .... 


16 


57 


58 


-45 


56 


36 





26 


9 


69 


7 


41 








7 


1 


NGC 6268 


17 


02 


40 


-39 


44 


18 





31 


10 


19 


8 


50 


1 





2 


22,34 


NGC 6322° . . . 


17 


18 


30 


-42 


56 


00 





50 


10 


50 


7 


06 


1 


8 


43 


1,18 


NGC 6425 .... 


17 


47 


02 


-31 


30 


00 





30 


9 


89 


7 


35 








25 


1,29 


NGC 6530° . . . 


18 


04 


48 


-24 


20 


00 





25 


10 


62 


6 


87 


1 


2 


> 41 


1 


NGC 6531 .... 


18 


04 


36 


-22 


30 


00 





20 


10 


47 


7 


07 


1 





49 


1,9 


NGC 6604° . . . 


18 


18 


06 


-12 


14 


00 





72 


11 


15 


6 


81 





7 


> 63 


1 


NGC 6613 .... 


18 


19 


54 


-17 


08 


00 





31 


10 


48 


7 


90 


1 





29 


12,14,34 


NGC 6664 .... 


18 


36 


42 


-08 


13 


00 





56 


10 


70 


7 


66 


2 


2 


36 


24 


Ruprecht 79 b . 


09 


40 


59 


-53 


51 


00 





53 


12 


55 


7 


65 


1 


31 


288 


1,28,32 


Ruprccht 119 . 


16 


28 


15 


-51 


30 


00 





42 


11 


64 


7 


18 


1 


9 


50 


1,22,34 


Ruprecht 127 b 


17 


37 


51 


-36 


18 


00 





77 


10 


92 


7 


35 





18 


97 


1,18 


Ruprecht 140 b 


18 


21 


51 


-33 


13 


00 





37 


12 


35 


7 


50 





1 


25 


34 


Stock 13 c 


11 


13 


05 


-58 


53 


00 





16 


11 


63 


7 


22 


1 


1 


64 


1,17 


Stock 14 


11 


43 


48 


-62 


31 


00 





19 


11 


66 


7 


06 


1 


1 


47 


1,14 


Trumplcr 7 . . . . 


07 


27 


22 


-23 


57 


00 





21 


11 


02 


7 


40 





5 


24 


14 


Trumpler 18 . . 


11 


11 


28 


-60 


40 


00 





23 


10 


66 


7 


19 





5 


71 


1 


Trumpler 20 . . 


12 


39 


34 


-60 


37 


00 





26 


11 


92 


8 


20 





1 


10 


34 


Trumpler 27 a . 


17 


36 


20 


-33 


31 


00 





89 


10 


42 


7 


06 


1 


2 


> 34 


1,29 


Trumpler 28 . . 


17 


37 


00 


-32 


29 


00 





45 


10 


88 


8 


30 








19 


14,34 


Trumplcr 34 . . 


18 


39 


48 


-08 


25 


00 





56 


11 


20 


8 


10 


1 


1 


18 


34 


vdB-Hagcn 217 


17 


16 


06 


-40 


49 


00 





60 


10 


42 


7 


54 








12 


2,34 



a Clustcr parameters unreliable. 
b May not be a true cluster. 
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c May contain pre-MS Herbig Be stars. 



References. — (1) WEBDA: (2) Ahumada et al. 2000; (3) Alfaro & Delgado 1991; (4) Baumc et al. 2003; (5) Carraro & Munari 2004; (6) 
Delgado, Alfaro, fc Cabrera-Cano 1997; (7) using mean values from FitzGcrald et al. 1979; (8) Forbes & Short 1996; (9) Hagen 1970: (10) 
Herbst & Havlen 1977; (11) Kalirai et al. 2001; (12) Lindoff 1971; (13) Lloyd Evans 1969; (14) Lynga 1987; (15) Moffat & Vogt 1973; (16) 
Moffat & Vogt 1975a; (17) Moffat & Vogt 1975b; (18) Moffat & Vogt 1975c; (19) Paunzen & Maitzen 2001; (20) Perry, Hill, & Christodoulou 
1991; (21) Piatti, Claria, & Bica 2000a; (22) Piatti, Bica, fc Claria 2000b; (23) Sanner et al. 2001 (24) Schmidt 1982; (25) Shobbrook 1985; 
(26) Subramaniam & Sagar 1999; (27) Stetson 1981; (28) Tadross 2001; (29) The & Stokes 1970; (30) Vazquez et al. 1997; (31) Vogt & Moffat 
1972; (32) Walker 1987; (33) Wesselink 1969; (34) this work. 
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Table 3 
Cluster mass distribution 



Cluster 


r 


N tot 


N ■ ■ 

iv missing 


AM 


a (M G ) 


Hogg 22 . . 


-1.11 ±0.27 


49 ± 10 


27 


3.9 


- 24.8 


NGC 3293 


-1.00 ±0.14 


189 ±21 


98 


4.2 


- 18.4 


NGC 3766 


-1.67 ±0.10 


191 ± 9 


45 


2.7 


- 12.2 


NGC 4755 


-1.04 ±0.13 


247 ± 17 


69 


2.7 


- 20.0 


NGC 6530 


-0.68 ±0.34 


66 ± 7 


25 


2.7 


- 13.8 


NGC 6604 


-1.27 ±0.19 


127 ±18 


64 


3.8 


- 20.8 



a AM represents the observed range of masses for the B stars in 
each cluster. 
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Table 4 

Distribution of open clusters 



Cluster 


I 

(deg) 


b 

(deg) 


d Q 
(kpc) 


dcent 

(kpc) 
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Table 4 — Continued 





I 


b 


d & 


dcent 


Cluster 


(deg) 


(deg) 


(kpc) 


(kpc) 


Stock 13 


290.52 


1.57 


2.12 


7.52 


Stock 14 


295.24 


-0.64 


2.15 


7.35 


Trumpler 7 . . . . 


238.28 


-3.38 


1.60 


8.95 


Trumpler 18 


290.98 


-0.14 


1.36 


7.62 


Trumpler 20 . . 


301.49 
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7.04 


Trumpler 28 . . 


355.99 
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1.50 


6.50 


Trumpler 34 . . 
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6.45 


vdB-Hagen 217 


346.76 
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1.21 


6.82 
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Table 5 

B Star Rotational Velocities 





M B 


Rb 


Vcrit 


(vsini) 




Spectral Type 


(M Q ) 


(Re) 


(km s- 1 ) 


(km s- 1 ) 


(usini) /v crit 


Bl 


13.5 


5.9 


560 


127 


0.23 


B4 


5.2 


3.2 


470 


108 


0.23 


B7 


3.5 


2.5 


440 


152 


0.35 


B9 


2.4 


2.1 


390 


134 


0.34 
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